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I briefly summarize the parallel sessions on Automated Calculation and Simulation Systems for high energy 
particle physics phenomenology at ACAT2002 (Moscow State University, June 2002) and present a short overview 
over the current status of the field and try to identify the important trends. 


1. Introduction 

Future particle colliders—planned or already 
under construction—will explore a new frontier 
in energy and precision. Final states with more 
tagged particles and better defined jets will be¬ 
come available for physics analysis. 

If there is low energy supersymmetry, the de¬ 
termination of the quantum numbers of the pre¬ 
dicted particles will require comprehensive stud¬ 
ies of complicated cascade decays, taking ad¬ 
vantage of spin correlations. If there is no 
low energy supersymmetry, physics beyond the 
standard model will contribute to processes like 
W/Z couplings in e+e“ —> / 1 / 2 / 3 / 4 , WW scat¬ 
tering in e+e" ^ 1 ^ 6 / 1 / 2 / 3/41 tt production 

in e+e“ ^ / 1 / 2 / 3 / 4 / 5/6 and WW/ti scattering 
in e+e“ ^ VeDefif 2 hfAf 5 fe, but the effects will 
be very small and can only be observed in preci¬ 
sion measurements of these processes—and com¬ 
plementary processes at hadron colliders—that 
must be compared to precision calculations in 
the standard model, effective field theories and 
specific models for physics beyond the standard 
model. 

Therefore, we will need reliable predictions, im¬ 
plemented in precision simulation tools to unleash 
the full potential of future colliders. Since the 
primary objective of the next generation of col¬ 
liders is the determination of the nature of elec- 
troweak symmetry breaking, particular care must 
be taken to obtain gauge invariant predictions. 
Polarization will be important to focus on longi¬ 
tudinal gauge bosons and for measuring quantum 
numbers of supersymmetric particles. 


It will be qualitatively more complicated to 
reach reliable high precision predictions for these 
multi particle final states than for the two par¬ 
ticle final states that have dominated physics at 
LEPl: 1) the parameter space of models beyond 
the standard model is much larger, 2) the num¬ 
ber of Feynman diagrams contributing to a pro¬ 
cess explodes combinatorially, 3) algebraic ex¬ 
pressions become much more complicated with 
the growing number of building blocks (i. e. in¬ 
dependent momenta and polarization vectors), 
4) gauge cancellations become numerically haz¬ 
ardous, and 5) phase space also becomes much 
more intricate. 

Therefore, even if we had enough graduate stu¬ 
dents and postdocs, we should not waste them on 
repetitive “assembly line” calculations. Instead, 
efforts must be directed toward formalizing the 
calculations so that the repetitive part can be 
delegated to computers, which are more patient 
with and particularly reliable for repetitive work. 
Of course, this is not a new observation and has 
formed the rationale for the AIHEP workshops 
starting in 1990, which were the precursor of the 
present AC AT conferences. 

In a bird’s eye’s view, the status and trends 
at ACAT2002 were: 1) tree-level calculations in 
the standard model for 2 ^ 4 and e“''e“ ^ 6 
have been well under control for some time now 
and can readily be used by non-experts, 2) the 
support for general 2^6 and 2 —*■ 8 processes 
is coming along, but not completely usable for 
production yet (systems with already complete 
physics support do not scale optimally, while opti¬ 
mally scaling systems have not been implemented 
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Figure 1. Components of the WHIZARD ACSS. 
WHIZARD can use different components for gen¬ 
erating matrix elements and relies on the VAMP 
library for adaptive multi channel sampling of 
phase space parameterizations matching the ma¬ 
trix element’s dominant peaking structures. The 
other ACSS are structured in a similar way. 


completely yet). 

The main challenges lie beyond the standard 
model (in the MSSM in particular) and vir¬ 
tual radiative corrections. Regarding the former, 
progress has been reported at the conference and 
regarding the latter, there is growing evidence, 
concern and consensus that the classic analyt¬ 
ical approach “does not scale”. Recently, sev¬ 
eral projects have started to work on different 
(semi-)numerical approaches and progress by one 
group has been reported at this conference Q. 


2. Automated Calculation and Simulation 
Systems 

Fully Automated Calculation and Simulation 
Systems (ACSS) in particle physics aim to pro¬ 
duce cross sections and event samples for a given 
physics model and scattering process 
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with as little as possible human intervention (ide¬ 
ally without any, of course). The job of ACSS can 
roughly be be divided into two steps 

1. calculate the matrix element T (i. e. gener¬ 
ate Feynman diagrams, derive arithmetical 
expression and generate executable code) 

2. integrate |Tp with a phase space measure 
d/r or generate events according to |T|^d/r 
(the efficiency of this step can be improved 
significantly by using information on the 
structure of T). 

In a third step, the partonic final states have to 
be transformed into hadronic final states. This 
fragmentation and hadronization is typically del¬ 
egated to the established players in this field like 
Pythia and HERWIG. 

Still, an ACSS is a complex piece of soft¬ 
ware and can only be constructed and maintained 
if broken into several components. As a typ¬ 
ical example, the structure of WHIZARD is dis¬ 
played in hgure |^. Indeed, some systems are 
complete and self-contained, while some provide 
components of complete systems]^: CompHEP is a 
complete system (two talks and two posters at 
ACAT2002), CalcHEP is a CompHEP clone (one 
talk), GRACE is a complete system (one talk 
and additional talks in the sessions on numeri¬ 
cal integration and sampling), O’Mega calculates 
efficient matrix elements for many particle fi¬ 
nal states, WHIZARD performs phase space inte¬ 
gration for many particle hnal states and inte¬ 
grates other components into a complete system, 
HELAC/PHEGAS calculates efficient standard model 
matrix elements and performs phase space in¬ 
tegration, Madgraph calculates standard model 

^Lack of space makes it impossible to provide complete 
references. The ordering proceeds from complete systems 
to components and is alphabetical in each group. 
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matrix elements^, Alpha calculates efficient stan¬ 
dard model matrix elements for many particle fi¬ 
nal states, FeynArts/FeynCalc/FormCalc calcu¬ 
lates standard model and MSSM loop diagrams, 
CalcPHEP (now SANG) calculates standard model 
2 —> 2 processes at one-loop (two talks and two 
posters). 

3. Trends at AC AT 2002 

3.1. Supporting More General Models 

Currently, all ACSS use a set of Feynman rules 
(i.e. a set of propagators and vertices) as in¬ 
put. This is mathematically equivalent to the 
Lagrangian describing the model. However, the 
derivation of the Feynman rules derivation from a 
Lagrangian—even if straightforward—can be ex¬ 
tremely tedious and error-prone (e. g. there are 
many thousand different vertices in the MSSM). 

In two talks tools for the derivation 

of Feynman rules in complicated models—the 
MSSM in particular—have been discussed at 
ACAT2002. Such tools need to be able to ex¬ 
pand lagrangians in terms of component fields 
and momenta and simultaneously handle mixing 
from non-diagonal mass matrices. These opera¬ 
tions must use a very natural notation in order to 
minimize the potential for human errors. 

Therefore, such tools evolve quickly into 
fully developed symbolic manipulation languages, 
with some special features: 1) transparent nota¬ 
tion for objects from mutually commutative, but 
non-commutative algebras, like non-simple gauge 
groups (this is poorly supported in general pur¬ 
pose computer algebra systems) and 2) support 
for the dual role played by covariant derivatives as 
operators and field components. Soon there will 
be several independently derived sets of machine 
readable Feynman rules for the MSSM. 

In addition, moving beyond the standard model 
requires that some systems have to relax assump¬ 
tions that were useful for implementing the Stan¬ 
dard Model in the first generation of ACSS: Ma- 
jorana spinors have to be supported in addi¬ 
tion to Dirac spinors, vertices of degree higher 
than four, more complicated momentum depen- 

companion program MadEvent |0] for event generation 
was introduced not too long after tne conference. 


dence in higher dimensional operators, non-local 
vertices factors (e. g. in non-commutative field 
theories), non Feynman diagram contributions 
(e. g. K-matrix unitarization), etc. 

One solution is to replace previously hard¬ 
coded subsystems making special assumptions by 
general purpose components. One poster re¬ 
ported how CompHEP’s special purpose C-routines 
for the symbolic evaluation of squared diagram 
are being replaced by calls to the general purpose 
form system. 

3.2. Component Architecture and Persis¬ 
tency 

The components in an ACSS have to commu¬ 
nicate across a variety of boundaries: abstraction 
layers (subroutines, objects, modules, etc.), ad¬ 
dress spaces (processes, networks, etc.), and time 
(data storage, persistency). These issues were ad¬ 
dressed by two talks . 

Pythia7 H is not an ACSS in the the strict 
sense, but all ACSS have to talk to PYTHIA or 
HERWIG for fragmentation and hadronization of 
hard events. PythiaT is a project to completely 
rewrite the Lund event generators in C++. It 
will provide a general structure for implementing 
models for event generation, not only the Lund 
Model (in fact, HERWIG++ has joined the effort). 
Today, PythiaY exists as a proof-of-concept ver¬ 
sion (some basic 2 —> 2 matrix elements, rem¬ 
nant handling and Lund string fragmentation. A 
new pre-release with the first HERWIG contribu¬ 
tions is planned for 2002 and a usable generator 
is planned for the following year. PythiaY should 
be available as the standard generator at the LHC 
start-up. 

External representations of event samples are 
required for communicating partonic events to 
soft QCD Monte Carlos and to detector simula¬ 
tions. External representations are also useful for 
saving both intermediate and final results. The 
latter can be useful if reading the stored events 
is faster than re-generation. However, there will 
always be a gain for reweighting processed data 
and one must also not underestimate the value of 
external representations for debugging. 

A proposal of a text based data description lan¬ 
guage for event samples and the implementation 
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of the accompanying toolkit was welcomed, 
but during the discussion there was demand for 
replacing the concrete syntax by XML, while re¬ 
taining the toolkit interface, was voiced. 

3.3. Loops 

In two talks 0,^ and two posters, the project 
of a re-implementation of the complete standard 
model one-loop radiative corrections (inch soft 
bremsstrahlung) for e“'"e“ —> // from scratch was 
presented. 

One motivation for such a re-implementation 
lies in the need to preserve the body of knowledge 
for future generations by providing a consistent 
and systematic option for redoing the calcula¬ 
tions. An application will be the relaxation of Z- 
pole assumptions for applications to off-resonance 
physics at a linear collider. 

The calculational procedure starts by a de¬ 
composition into form factors, followed by the 
calculation of the form factors using a collec¬ 
tion of forms procedures for symbolic calcula¬ 
tion, renormalization and generating numerical 
Fortran code. Currently, the t-channel neutral 
current processes and and the s-channel processes 
are done and show very good agreement with ex¬ 
isting codes. Interference contributions and de¬ 
cays are in progress. 

4. Outlook 
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Figure 2. Schematic representation of the sta¬ 
tus of fully automated calculation tools. (The 
names are a random representative of a class of 
programs.) Important contributions can be made 
in the “terra incognita”, that is covered neither by 
the existing many-loops/few-legs nor the many- 
legs/no-loops programs. 


among developers and among programs—using 
pluggable components in the ideal case—is nec¬ 
essary for this. 
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If we would ever manage to produce a fully gen¬ 
eral and efficient ACSS, that is even easy to use 
for non-experts, we would risk to make ourselves 
obsolete as phenomenological theorists. Fortu¬ 
nately, we will never get there because theoretical 
theorists will always come up with new features 
that we have not anticipated and experimentalists 
will always push the frontiers in energy and pre¬ 
cision, calling for ever more precise calculations 
of ever more complicated processes. 

Currently there is a vast uncharted territory 
(figure ^ and a systematic exploration of the high 
energy/high precision region with many legs and 
many loops will require the development of new 
semi-numerical methods. 

Finally, variety is good, since we need to be 
able to cross check our results. Communication 
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